Stability, complex modes and non-separability in rotating quadratic potentials 
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We examine the dynamics of a particle in a general rotating quadratic potential, not necessarily 
stable or isotropic, using a general complex mode formalism. The problem is equivalent to that of a 
charged particle in a quadratic potential in the presence of a uniform magnetic field. It is shown that 
the unstable system exhibits a rich structure, with complex normal modes as well as non-standard 
modes of evolution characterized by equations of motion which cannot be decoupled (non-separable 
cases). It is also shown that in some unstable cases the dynamics can be stabilized by increasing 
the magnetic field or tuning the rotational frequency, giving rise to dynamical stability or instability 
windows. The evolution in general non-diagonalizable cases is as well discussed. 

PACS numbers: 03.65.Fd,03.65.Ca,03.75.Kk 



INTRODUCTION 



Quadratic forms in boson operators or generalized co- 
ordinates and momenta are an ubiquitous presence in the 
theoretical description of diverse physical systems. They 
often arise through the linearization of the equations of 
motion around a stationary point, as in the case of the 
random phase approximation (RPA) [l], Q , providing a 
basic tractable scenario. They play an important role 
in the description of Bose-Einstein condensates (BEC) 
as well as in other areas like quantum optics PI, 
disordered systems [l(| and dynamical systems pll - [l3| . 
Nonetheless, while positive forms, characteristic of sta- 
ble systems, are well known to be diagonalizable, i.e., 
they can be written in terms of normal coordinates and 
viewed as a set of independent bosons or separate oscil- 
lators (plus eventually free particle terms in the presence 
of standard zero frequency modes 0, 0]), non-positive 
ones may not admit such a diagonal representation [l4j ]. 
Non-positive forms can arise in the description of BEC 
instabilities and fast rotating condensates [l5l - l20| . 

as well as in generalized RPA treatments [2l], [22| • 

In [3] we have extended the standard methodology 
for diagonalizing quadratic bosonic forms by using gen- 
eralized quasiparticle operators fulfilling boson-like com- 
mutation relations, associated with non-hermitian coor- 
dinates and momenta. This allows to characterize the 
operators exhibiting purely exponential evolutions (com- 
plex modes) in non-positive forms, enabling a precise de- 
scription of the dynamics and quadratic invariants. We 
have also pointed out that non-positive forms can in 
some cases be dynamically stable, as the evolution can 
remain quasiperiodic, irrespective of the initial condi- 
tions. Moreover, we have noticed the existence of non- 
diagonalizable cases where the equations of motion can- 
not be fully decoupled and which may arise even if all 
eigenfrequencies are non-zero. The method of ref. [f4j 
has been found useful in the context of BEC, being em- 
ployed to study the emergence of instabilities in trapped 
BEC with a highly quantized vortex through the 

Bogoliubov-de Gennes equations. 

Here we will apply this methodology to the basic prob- 



lem of a particle in a rotating anisotropic quadratic po- 
tential, not necessarily stable. This system is formally 
equivalent to that of a charged particle in a uniform mag- 
netic field in a general quadratic potential [HI, [HI, [23| . 
The problem is therefore relevant for many fields. In 
particular, the rotating case has recently become rele- 
vant in the context of BEC in rotating anisotropic traps 
[T5M20I ] . which in the Landau level approach are basi- 
cally described by a cranked quadratic potential of the 
type here considered. The stable system is well known 
[3, H, [H, [24| and diverse aspects of the stable anisotropic 
rotating case in the context of rotating condensates have 
recently been investigated [f7l - [2fl |. 



We will here examine the general unstable case, which 
is of interest for fast rotating condensates as the Hamil- 
tonian ceases to be positive definite at high frequency 
due to the centrifugal force. As we shall see, the unsta- 
ble system exhibits a rich structure, with several differ- 
ent dynamical regimes as well as some quite remarkable 
features, including: a) The possibility of becoming non- 
separable at the boundaries of regions with distinct dy- 
namics, in the sense that the Hamiltonian can no longer 
be written as a sum of two independent standard or com- 
plex modes. In such cases the system will exhibit anoma- 
lous evolutions characterized by a set of linear equations 
which cannot be decoupled and which may lead to coor- 
dinates and/or momenta evolving with terms cx te xt or 
even oc t 3 ; b) The possibility of achieving dynamical sta- 
bility in some unstable cases by increasing the magnetic 
field or tuning the rotational frequency. In particular, 
a stable anisotropic rotating potential becomes dynam- 
ically unstable just in a finite frequency window, recov- 
ering dynamical stability at high rotational frequency, 
whereas an unstable saddle type potential can become 
dynamically stable in a certain frequency window. 



In sec. II we briefly revisit the main features of the 
formalism, discussing the concept of separability in gen- 
eralized coordinates and momenta and the evolution for 
general non-separable cases. The application is discussed 
in III while conclusions are drawn in IV. 
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II. FORMALISM 



A. General evolution and dynamical stability 



We consider a general quadratic hamiltonian 

h = | TjjPiPj + Vjjgjgj + UjjjqiPj + Pjgj) (la) 



\R t T-L c R , i? — ( ^ ) , H c — f j, 



(lb) 



where i denotes transpose, T, V are symmetric matrices 
and p, gr are hermitian coordinates and momenta satisfy- 
ing \Pi,Qj] = -iSij, [qi,q 3 ] = \pi,Pj] = 0, i.e., 



RR l - (RR'y = M c , M c = i I _j 



1 




(2) 



The ensuing Hcisenberg equations of motion lead to a 
closed set of linear equations which can be written as 



= -[h,R]=U c R, 



dt 

H c = MM = i 



U* T 

-V -U 



(3) 
(4) 



where T-L c represents the RPA matrix in coordinate rep- 
resentation [14]]. It completely determines the system 
dynamics. Its eigenvalues come in pairs of opposite sign 
and can be complex in unstable systems. Moreover, R c 
can also be non- diagonalizable (as in the case of free par- 
ticles U = V = 0, Tij — tiSij, although other cases can 
also arise, as discussed later). A positive definite R c 
(Rt^icR > V real R ^ 0) ensures a diagonalizable T-L c 
together with a real spectrum [J, 0, EH (standard stable 
case) but the converse is not true. 

Under a general linear canonical transformation R = 
UR' , with the matrix U satisfying UMcU* — M c in 
order to preserve Eq. @, we have H' c = U t 'H c L( but 
H' c = MM' C = U^RJA, ensuring the invariance of the 
eigenvalues and the Jordan canonical form of R c . Ma- 
trices O = M c O are precisely those accounting for the 
closed algebra of the forms ((J) : If Oi 



\R t O i R, then 



[Oi,Oj] = \R*CR, with C = [O^Oj] 



(5) 



It is obviously equivalent to use a representation of h in 
terms of boson operators bj ,b'- — (qj ± ipj ) / y/2 satisfying 
= Sij- Defining Z = (£ t ) and the unitary matrix 
S = ;4| (i -i), such that Z = SR, we may rewrite h as 



i,3 



\Z^RZ, H — SH C S^ — ( 



For a time independent h, the solution of system ([3]) is 

R(t) = exp[-iH c t]R, (7) 

where R = R(0). Eq. ([7]) is itself a linear canonical trans- 
formation. A system which is dynamically stable, i.e., 
leading to a bounded quasiperiodic evolution of all op- 
erators Pi, qi, corresponds to a matrix H c which i) is 
diagonalizable and ii) has only real eigenvalues: 

a) If He is diagonalizable, such that R c = W^W -1 
with {R! c )^ v = A„<5 M „, we may expand Eq. (J7J as 



(6) 



A = ~\V + T -i(U -U% B ± = ±[V -T±i(U + U*)] 

The ensuing RPA matrix, defined by idZ/dt = RZ, is 
just R = SH C S> and has obviously the same eigenvalues 
(and Jordan canonical form) as R c . 



R(t) 



-iX u t 



(8) 



where W„ is the v th column of the eigenvector matrix VV 
and Z' = W _1 i? is a set of normal operators satisfying 

rl7' 

"=KZ' V , (9) 



-%\ v t 



Z' v . For complex 



dt 

and evolving then as Z' u (t) 
eigenvalues A„ , these generalized normal operators repre- 
sent exponentially increasing or decreasing modes (com- 
plex modes) and the dynamics is unbounded (±A„ are 
both eigenvalues). They can always be ordered and nor- 
malized such that those associated with ±A„ (b' v and b' v , 
with Z' = (b' ,b') r ) satisfy boson-like commutation rela- 
tions 0, i.e., = 5^, [b'^X] = [V„,V V ] = 0, but 

b' v 7^ b'l if A„ is non-real (l4|. 

b) If R c is non-diagonalizable, the system @ cannot 
be fully decoupled, but we may use its Jordan canonical 
form H c = WR^W" 1 , with H' c having blocks of the form 



l~i! r — 



A„ 1 
\ v 1 



V 






K 



\ 



We may then expand Eq. |7J) as 



R(t) 



^E 
fe=i 



E z '"i 

l=k 



d-k 



(I k)\ 



(10) 



where again Z' 
and i>k , k = 1 , . 



W d v is the dimension of the block 
. , d v labels elements within each block, 



with Wi/ k the Vk column of the generalized eigenvector 
matrix W. The generalized normal operators Z' satisfy 
the "minimally coupled" evolution equations allowed by 
the Jordan form, i.e., idZ'/dt = R' C Z' or 

dZ' 

5k, d „)Z', (11) 



dt 



(1 



leading to Z' Vk (t) = EZk ^^T- 



The dynamics 

is then unbounded even for real \ v . The free particle 
case corresponds to \ v — and d v — 2. Other cases are 
discussed in sec. III. 
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B. Separability 

We will denote the Hamiltonian ([1]) as separable if 
there is a linear canonical transformation R = UR! (with 
UMJA 1 = M c ) such that H' c = U l H c U is diagonal. In 
this case we may then rewrite h as a sum of independent 
elementary quadratic systems, 

h = ±]Tk^ + /W*), (12) 

where [g^pJJ = i^, K,^] = [p^,P^] = 0. The 
diagonal form (fT2|) is not unique, as p' u , q' v can be 
rescaled ((p' u ,q' u ) -> (ap' u ,q' u /a)) or swapped ((p' v ,q' v ) -> 
(ill—pi)), but the products a v p v — \ 2 determine the 
eigenvalues of T-L c and are hence unique. 

In contrast with the conventional normal mode expan- 
sion of a positive definite h, each of the terms in (|12[) 
can here represent not only i) a standard stable oscil- 
lator (a v > 0, fi u > 0), but also ii) an "inverted" os- 
cillator (a„ < 0, fi v < 0), iii) a generalized free particle 
(a„/3„ = 0, with a v ^ or j3 v ^ 0, which can be standard 
or inverted), iv) an "unstable oscillator" (a v j3 v < 0) and 
v) a "complex oscillator" (a. v @ v complex), where p' v , q' v 
are no longer hermitian (IA complex). We should also add 
the vanishing case (0) a v = j3 u = 0, where boihp' v and q' v 
commute with h. Separability in hermitian coordinates 
and momenta (U real) is a restricted class of separability 
since the eigenvalues ±y/a v [} v of i-L c are in such a 
case real or imaginary, while diagonalizable cases with 
full complex A„ do exist [l4| . 

A diagonalizable H c ensures separability since in this 
case we may rewrite h in terms of the generalized normal 
operators Z' = (&', 6')* of Eq. (jgj) as fl3| 

h = ^A,(&X+|) = IE A ^ + <?'*), ( 13 ) 

where q' v = {b' v + V v )/y/2, p' v = (b' v - V v )/(y/2i), i.e., 
U = WS. For complex A„, b' v ^ b'\ and p'„,ql are non- 
hermitian. In case iv) we may still rewrite the ensuing 
term in (|13[) in terms of hermitian q' u , p' v by a complex 
rescaling {p' v ,q' v ) — > (pl/Vi,Viql) [IJ], while in ii) we 
should choose A„ < for p' v , q' v hermitian. Conversely, 
for a v (3 v ^ or a v — (3 U — 0, each term in (fT2"j) leads to 
a diagonalizable 2x2 block in H' c . However, the separa- 
ble case also includes the free particle case iii) where T-L c 
is non-diagonalizable, as T-l' c will contain a Jordan block 
Ti' Q = (q q) with \ u = 0. Here H 2 remains diagonalizable. 

Hence, systems where W 2 is non-diagonalizable, im- 
plying H! c having a Jordan block of dimension d v > 2 or 
d v = 2 and A^ =^ 0, are non-separable. They may arise 
even in simple unstable cases (sec. Ill) and their evolution 
can be determined through the general solution (fTOj) . 

Dynamically stable quadratic systems correspond to a 
separable h with terms just of the form i), ii) or 0), and 
have then a discrete spectrum 

^{n„} = E,A,K + |), (14) 



where A w > 0, < or in cases i), ii) or 0). An example 
of a non-positive dynamically stable form is an angular 
momentum component / = q+P- — q~p+ = \ (p 1 + +q' + ) — 
Up'- + where p'± = p± - q T /2, q' ± = q±/2+p' T , 

which is the sum of a positive plus an inverted oscillator. 
Here C c is diagonalizable with eigenvalues ±1, 

III. APPLICATION 

We will consider the quantum problem of a particle in 
an anisotropic quadratic potential, not necessarily sta- 
ble, rotating around one of its principal axes (z). It 
is formally equivalent to that of a particle of charge e 
in a uniform magnetic field H parallel to this axis in a 
quadratic potential. The Hamiltonian of the latter reads 

H = (P ~* A/C)2 + \{K X X 2 + K y Y " + ^ 2 )(15) 
2m 

p2 + p2 

= \ [ — v + K' X X 2 + K' y Y 2 - nL z ] + H z (16) 

where A — \H x R is the vector potential, il = e\H\/mc 
the cyclotron frequency, L z — XP y — YP X the angular 

momentum component, H z — \{— + K Z Z 2 ) and 
K, y = K x , y + mfi 2 /4 . 

For fi — > 20, Eq. (|T6|) is just the cranked Hamiltonian 
describing the intrinsic motion of a particle in a rotat- 
ing quadratic potential with constants K' x y : If H(t) = 
U(t)H(0)W(t), with U(t) = e -^LA/n^ t h e Heisenberg 
equations for rotating operators 0(t) = U(t)OU' (t) are 
those for the ^-independent cranked Hamiltonian H = 
H(0) - VLL Z . 

Since H z is fully decoupled from the rest and its treat- 
ment is trivial, it will be omitted in what follows and 
all considerations will refer to the motion in the x, y 
plane. Defining dimensionless operators q — R^J mflo/h, 
p = P/^hmflo satisfying [g^jy] = i5^', where Qq is a 
reference frequency, we can rewrite h = H xv /MIq as 

h = \[pl+pl + k' x q 2 x + k y q 2 y ]-ul z (17) 
= ^R'HcR, k x ^ y = k x , y +uj 2 , (18) 

with fc M = ^/(mfig), uj = f2/(2f2 ), h = q x p y - q y Px 
and i?* = {q x ,q v ,p x ,p y ). 

From the form of Eq. (|15[) , it is apparent that for fixed 
k^, the field cannot change the number of positive or 
negative eigenvalues of the Hamiltonian matrix % c (the 
number of positive and negative diagonal elements is the 
same in any real diagonal representation of a quadratic 
form). For m > 0, % c will then have none, one or two 
negative eigenvalues if and only if none, one or both of the 
constants k x ,k y are, respectively, negative. The positive 
definite case corresponds then to k x > 0, k y > V oj, i.e., 
k'^ > uj 2 0i 13 1 although we will now see that dynamical 
stability is not restricted to this case. 
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A. Dynamical stability 



The RPA matrix H c = M C H becomes 



1 





lo 


1 





— LU 








1 




k x 








LU 


\ 





K y 


— LU 






(19) 



and its eigenvalues, which are the system eigenfrequen- 
cies (in units of M7q), are ±A + ,±A_ with 



A+ = 



^(k' x + k> y )/2- 



±A, 



k' y ) 2 /A + 2u 2 (k' x +k' y ). 



(20) 
(21) 



They satisfy A+A?. = clctf'Hc] = k x k y . We first note that 
for k x k y ^ and A ^ 0, H c is diagonalizable since it will 
have four different eigenvalues. We will then show that 
for A = 0, H c is non- diagonalizable and, moreover, h is 



non-separable. At fixed kn, A 



if k XyV < and 



= = \\yf^k x ± \f~^k~ v 



(22) 



It also follows that A± can be both real only if k x k y > 0. 
For a charged particle in a magnetic field, this opens 
up the possibility of full dynamical stability around a 
quadratic maximum (k Xiy < 0) but dismisses it for a 
saddle point (k x k y < 0). It is indeed verified that A± are 
both real for fc x „ > as well as for k x _„ < and 



*x,y > as well as for k XtV 

M > ^ 



(23) 



At fixed k^, the dynamics in the vicinity of a quadratic 
maximum can then be stabilized by increasing the field. 
The behavior of A± for increasing \uj\ at fixed k^ is de- 
picted on the left panels of Fig. [1] As the scaled cyclotron 
frequency 2\lu\ is increased, at an anisotropic maximum 
A± evolve from imaginary (|w| < lu~) to full complex 
(lo~ < \lu\ < cj+) and finally to real (\lu\ > ljj) values, 
reaching the system dynamical stability for \lu\ > lu+, 
whereas for k x k y < 0, A+ is real but A_ is imaginary Vw. 

Stability in rotating potential. At fixed k'^, the previous 
picture is seen quite differently and leads to dynamical 
stability and instability windows in the anisotropic case 
k'x k' y (right panels in Fig. [1]), i.e., when the rotation 
has a non-trivial effect. Owing to the centrifugal force, 
H c is here positive definite just for lu 2 < k'. However, 
A± are real also for lo 2 > k' x y > 0, implying that motion 
in a rotating stable potential becomes dynamically stable 
at high rotational frequencies, and dynamically unstable 
just in the finite interval 



Mm[k' x ,k' y ] < lu 2 < Max[k' x , k' y ] 



(24) 



where A_ becomes imaginary or (see III B). In contrast, 
for a saddle point with k' x > > k' y (or viceversa) the 
system becomes dynamically stable in the windows 

(-K<k' y <o), ^ (25) 

(—3k' x < k' y < —k x ) , 



K <^ 2 

k' x <ijj 2 < lu' 2 



3 



1 1 1 — 

- k x < 


f 


■k' v /k'„ = l 5 

y / x 

: k' x >o 




_ 


■ sr. > 


i 








: k y /k x =-0.5 
- k x > 


i 1 y • 


■k'„/k'v=-l 5 
1 k' x >0 










\ 






i i i i i i i 


i 1 i 





FIG. 1: The real (solid, blue lines) and imaginary (dashed, 
red lines) parts of the eigenfrequencies (|20[) for selected fixed 
values of k XtV (left panels, corresponding to a particle in a 
magnetic field) and k' x>y (right panels, corresponding to a par- 
ticle in a rotating potential), in terms of the scaled cyclotron 
or rotational frequency. Vertical dotted lines separate dynam- 
ically stable and unstable regions. We have set cjo = V^PM 
(yklc) in the left (right) panels. 



with tut = 



\K-K\ 



the single value of \lu\ where A = 

at fixed k'^ (requires k' x + k' y < 0). Thus, l z can turn 
unstable the dynamics of a stable anisotropic oscillator 
and stabilize that around a saddle point within the limits 
(|25p. It can never stabilize a quadratic maximum. 
In the isotropic case k' x = k' y = k' , [h, l z ] — and Eq. 

(l20l) leads then to |A±| = \Vk? ±lu\. The rotation has 
no effect except for the shift ±lu. The instability window 
collapses into a single point lu 2 = k' where the dynamics 
remains stable (see below), since it corresponds to the 
Landau case k x = k y = [23|, [2f| . 



B. Separability 

Let us now examine the separable representation of 
Eq. (fT7|) . feasible for A ^ 0, and the ensuing distinct 
dynamical regimes. Defining, as in the stable case [l|,l2j, 



P± 



where 7 = 



q± = 



VPy,x 



777 



(26) 



2A-k'+k' 



— 27 



and [q v ,Pp] = i5 p 



4 W ) •! — k'^+k'y VHV^Hi — ""tiv, 

\p»,Pv] = [Qn,qu] = for [i,v = ± (with 7,77 ->• if 
lo — > and k' x > k' y ), we may rewrite (|T7)) for A ^ as a 
sum of two independent elementary quadratic forms , 

h - \{a +P 2 + + (3+ql) + \{a_p 2 _ + f3.q 2 _) , (27) 

- ^-^^ ,(28) 



a± = 



2A + k' !c ~k' v ±4uj 2 
4A 



with a±(3± = Aj_, For real A 7^ 0, p±, q± are hermi- 
tian and X± is real or imaginary. Eqs. (I2~?rf2"5)) are, how- 
ever, also applicable for imaginary A / 0, where p±,q± 



5 




4 - 



3 „ " 



_ 1 1 1 1 1 1 1 1 1 

D 

i t 


| i i | i i i | i _ 

! A 

.If ; 


E 

1 k 
■ 1 i i i 1 i i i 


~^\h a 

I D 

i , , i , , , i , - 



k' x /0J 2 

FIG. 2: Top: Regions in the k x ,y plane at fixed to with distinct 
dynamical regimes: A is the positive definite sector, F the 
Landau point, B the non-positive dynamically stable sector, 
C and D are unstable regions with four and two imaginary 
eigenfrequecies, E that with full complex eigenfrequencies, 
while the solid lines j, k depict the non- separability curve A = 
0, with degenerate real (j) or imaginary (k) eigenfrequencies 
(see text). The points L (k x = 0, k v = —4a; 2 or vice- versa) 
are the exceptional non-separable cases with a single vanishing 
eigenfrequency. Dashed lines indicate separable cases with a 
free particle term (standard in g,i and inverted in h) and a 
stable (g,h) or unstable (i) oscillator. Bottom: Same details 
in the k' Xty plane, corresponding to the rotating system: All 
regions and curves are just shifted by lj 2 . Plots in Fig. [T] 
depict the behavior of X± along straight lines running from 
co (cj = 0) to the origin (|uj — > oo). 



are non-hermitian and X± complex. From the previously 
mentioned property of diagonal quadratic forms, it fol- 
lows that all four coefficients a±, (3± must be positive 
for k x ^ y > 0, just one (two) of them will be negative for 
k x k y < (k x ,y < and A real), one of them will vanish 
for k y = 0, k x ^ (or vice- versa) whereas two of them 
will vanish in the Landau case k x = k y = 0. 

Diagonalizable cases. For A^0 and k x k y 0, X± =/= 
and ii c is diagonalizable. Defining p' ± — \/a±/X±p±, 
q' ± = vSM±9±i we ma y rewrite ([27]) as 

2 ,)=E A ^'X + 5)> ( 29 ) 



where b[ 



~ V2 



_ 1 



-ip 



V2 



are the generalized 



normal operators evolving as b' v (t) — e~ lXvt b' u , b' u (t) — 
can be directly obtained from the eigen- 
At fixed k^, the diagonalizable sectors are 



C Oj,, Willi 

vectors of H 



(fig.©: 

A) k XlV > 0: Here a± > 0, (3± > 0, with A± > 0. This is 
the positive definite case (case i in sec. II B). 

B) k x ^ y < 0, \uj\ > Here a + > 0, (3 + > but a_ < 0, 
(3- < 0, implying A + > but A_ < for p'_,q' _ hcr- 
mitian. Eq. (|29|) becomes a standard plus an inverted 
oscillator (cases i+ii), remaining dynamically stable. 

C) k XtV < 0, \to\ < bj~: Here a„/3„ < for v — ± and 
A± are both imaginary. Both terms in (|27[) are unstable 
oscillators (case iv), leading to p' ± , q' ± non-hermitian. 

D) k x k y < 0: Here a + > 0,/3+ > but a_/3_ < 0, with 
A+ real, A_ imaginary. Eq. p^|) becomes a stable plus an 
unstable oscillator (i+iv), withp'_,g'_ non-hermitian. 
In all previous cases p± , q± in (1271) are hermitian. 

E) k XyV < 0, uj~ < \w\ < o;+: Here a±, (3± and A± are full 
complex and p± , q± as well as p' ± , q' ± are non-hermitian 
(case v). They represent complex normal modes. 

F) k x = k y = (Landau case): Here a+ — 1, /3+ = 4cj 2 
whereas a_ = j3- = 0, leading to A + = 2\ui\ and A_ = 0. 
h is then a standard plus a vanishing oscillator (cases i + 
0). This well known case [23|,[25[ is then dynamically sta- 
ble in the x, y plane despite the vanishing egenfrequency. 

In cases A, B and F, p' u , q' v are hermitian, with b' v = 
b'\, and h possesses then a discrete spectrum 



X+n + + A_n_ + i(A+ + A_) (30) 



with A+ > 0, while A_ > in A, A_ < in B and A_ = 
in F. These are the dynamically stable cases. 

At fixed k'^ (rotating potential) all regions are just 
shifted by +us 2 (lower panel in fig. This shift leads to 
the different behavior of A± with u> depicted in fig. Q] 
Separable non- diagonalizable cases. They arise for 
^ or viceversa (sectors g,h,i in Fig. 
A = 2cj 2 + k x /2 and we obtain A_ = 
^ for A ^ 0, with 



k y = and k a 
For k y = 0. 



but A+ 



h = W. 



+ \W + 



i &x 9 

2 Al P - 



(31) 



In g, k x > and Eq. (f3"Tj) corresponds to a stable oscil- 
lator (A + > 0) plus a free particle (cases i+iii). In h, 
— AliS 1 < k x < and A+ is still real, but the second term 
in (|31[) becomes negative: h becomes a stable oscillator 
plus an "inverted" free particle term. The latter "ab- 
sorbs" here the instability, allowing dynamical stability in 
the coordinates p +1 q + . Finally, in sector i (k x < — 4u> 2 ) 
A+ becomes imaginary: Eq. pip corresponds here to an 
unstable oscillator plus a standard free particle (iv+iii). 

Eq. dHJ) leads to an evolution b' + (t) = e iA +*5' + , b' + (t) = 
e~ iA +*6' + ,p_(t) =p-,q-(t) = (k x /X\)p-t + q-, which is 
characteristic of a non-diagonalizable % c with a canonical 
form (A ~ B denotes A = W'^BW) 



Ho 



( A+ 
-A+ 
1 

V o ooo 



(32) 
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C. Non-separability 

We now examine the very peculiar case A = and 
lj 7^ 0, where the eigenfrequencies become degenerate, 



A+ = 



(3K + k')(3k'+k'J 



= X, and both H c and H 2 . are 



non-diagonalizable. It occurs at the threshold for full 
complex solutions, i.e., |w| = uif (ui' c ) at fixed k^ {k') 
and corresponds at fixed w to a parabola rotated 7r/4 
with respect to the k XtV axes, with vertex at k xy = —uj 2 , 



k~ 



(curves j, k, and point L in fig. [2]). 



j) Here \ui\ = w+ and A = ^/k x k y = *J cu^ 2 — u c 2 > 0. 
He can be reduced to two non-trivial Jordan blocks 



He 



A 1 
A 
-A 1 
-A 



(33) 



which indicates non-separability. The transformation 



P± = 



q± = 



\f~k~x 



-_(p Xt y + U) C q ViX ) 



x,y U7 C 



+ A 2 



2A 2 



\1x,y + ,2 



A 2 



allows to express h in this case as 



h = |(p+ +P-) - Kl+P- - Q-P+) , 



(34) 

Py.x) (35) 

(36) 



which is the cranked hamiltonian for a rotating free par- 
ticle and is a basic non-separable form: The equations 
of motion dp±/dt = ±Ap T , dq±/dt — ±Ag T + p±, can- 
not be fully decoupled even though there is no vanishing 
eigenfrequency. More explicitly, defining the operators 
b'± = b' ± = ±^±+^ satisfying [b' v ,b'»] = <W 

Wvi b'fj] = [b' v , &JJ = 0, we can rewrite Eq. (|36]l as 

h = X(b' + b' + -b'_b'_)-b' + b'_. (37) 



.db'. 



The ensuing equations of motion, i-.. 
TA6j_ + 61-, correspond exactly to the Jordan form (|3"3")) 
and possess the general solution 



b' ± (t) = e* iXt b' ± , b' ± (t) = e ±lM (b' ± - itb' T ) 

which gives rise to unbounded spiral-like trajectories in 
the variables p±, q±. This is apparent from Eq. ([36[) : 
Since [I, h] = 0, where I — q+p- — q~p+, the evolution 
operator e~ lht — e *^* e -*('»+^0 i represents a rotation of 
frequency A applied to a free motion. Such motion cannot 
arise from a separable h: Although h is the sum of two 
commuting quadratic forms and the equations of motion 
can be reduced to two separate blocks, the operators in 
each block do not commute with those of the other. 



k) Here \oj\ = uj c and A = i\jk x k y 



UJ, 



is imaginary. This case arises for k x ^ k y and leads 



to the same canonical form (|33|) . Replacing w+ <H- uj c 
Eqs. ([3"3 ]) -([3"5 ]) leads here to 



p 2 .) - \X\{q+P- + q-P+) 



(38) 



where p± , q± are hermitian. With a complex scaling 
p_ — » ip_, q_ — > —iq- Eq. ([35]) becomes identical to 
Eq. (|36l) with an imaginary A. The remaining equa- 
tions remain then unchanged but lead to exponentially 
increasing or decreasing evolutions for the operators b' ± , 
b'± , corresponding to a "boost" applied to a free motion. 

L) This is an exceptional critical point where cases 
B,C,D,E,h,i,j and k merge. Here k y = 0, k x = —4lu 2 
(i.e., k' y = uj 2 , k' x = — 3ui 2 ) or vice- versa, implying 
\oj\ = cj+ = uj~ and A± =0: All four eigenvalues of 
H c vanish. Nonetheless, H c is still of rank 3, implying 
that it becomes similar to a full Jordan block (d v = 4): 



/0 1 0\ 

10 

1 

\ / 



(39) 



instead of two free particle blocks, as would appear from 
Eq. (|33[) for A — > 0. In this case the equations of motion 
cannot be even partially decoupled. The transformation 

P± = 8±(p Xj y+uq yiX ), q± = ^(3q XtV - u^py^O) 

with 5+ — 1, 5_ = 2, allows to rewrite h at this point as 

h = \p\ - uq + p- , (41) 



which is again a basic non-separable form: The ensuing 
equations of motion, = 0, = -Ji = p + , 

~tt~ = —uq+i exhibit the structure of the Jordan form 
(|39| for ui 7^ and lead to a polynomial evolution of third 
degree in t for g_ : 

P-(t) = P-, p+(t) =p++utp^, 

q+(t) = q + +p+t+±ujt 2 p-, . (42) 

q~(t) = q- — ujtq + — ^ujt 2 p + — ^uj 2 t 3 p- 

Coordinates q + , g_ experience then a constant and lin- 
early increasing acceleration respectively. In terms of the 
operators b'_ = —ujp-, b' + = ip + , b' + = q+, b'_ = iq-/ui, 
which satisfy [b^b 1 ^] = S^ v , [b'^b v 
may also express Eq. (|4"Tj) as 



[b'^b'u] = 0, we 



b'b' 



hb'i 



(43) 



Their equations of motion, db'_/dt = 0, idb' + /dt = b'_, 
idb'+/dt = b'_i_, idb'_/dt = b' + , follow exactly the Jordan 
form (j3"9"| . Note that the disappearance of one of the ki- 
netic terms in (|41[) is not exceptional for a non-separable 
form: Eq. ([36]) (and hence ([38]) ) can also be rewritten as 
p\ — XI if q± —> q± + -r^p T . We may also rewrite (j4T|) 
with two kinetic terms with a similar transformation. 

Eq. ([391) suggests that this case could be considered as 
a "free inseparable pair" , generalizing the free particle 
case where H c ~ (o o)- For a free particle H 2 = while 
here Ht. = 0. 
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IV. CONCLUSIONS 

We have first analyzed, within the formalism of ref. 
[lij , the dynamics in general unstable quadratic bosonic 
forms, discussing the treatment of the general non- 
diagonalizable case and determining the conditions for 
dynamical stability and separability. We have then ap- 
plied the formalism to the basic problem of a particle 
in a general rotating quadratic potential, relevant in the 
context of fast rotating condensates in harmonic traps 
and formally equivalent to that of a charged particle in 
a uniform magnetic field in a quadratic potential. The 
present analysis unveils the rich variety of behaviors that 
can be exhibited by the unstable system, summarized 
in fig. [21 together with some quite remarkable features, 
which could lead to observable effects in fast rotating 
condensates. In particular, we have determined: a) The 
regions of dynamical stability. Intrinsic motion in a rotat- 
ing stable potential remains dynamically stable at high 
frequencies, becoming unstable just in a finite frequency 
window in the anisotropic case, whereas in a rotating 



saddle potential it can become dynamically stable in a 
certain window (Eqs. ([211), and fig-©; b) The re- 
gions in parameter space where H xy can written as a sum 
of two independent quadratic systems (separability) , em- 
ploying non-hermitian normal coordinates and momenta 
if necessary (sector E), and those where such a represen- 
tation is not feasible (non-separability); c) The explicit 
transformations and final forms for all cases, including 
the energy spectrum in the dynamically stable cases and 
the "minimally coupled" standard forms and equations 
of motion in the non-separable cases; d) The existence of 
an exceptional non-separable zero mode case (point L in 
fig. [2j) where all eigenfrequencies vanish and the Jordan 
Block has dimension 4. It is not equivalent to a standard 
zero frequency mode and leads to coordinates evolving as 
a third degree polynomial in time. These results indicate 
that similar peculiar effects can arise in more complex 
unstable quadratic systems, which can be analyzed with 
the same general formalism and techniques of sec. II. 

RR and AMK are supported by CIC of Argentina. 
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